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The TRAPPIST-1 system is remarkable for its seven plan- 
ets that are similar in size, mass, density, and stellar heating 
to the rocky planets Venus, Earth, and Mars in our own Solar 
System [1]. All TRAPPIST-1 planets have been observed with 
the transmission spectroscopy technique using the Hubble or 
Spitzer Space Telescopes, but no atmospheric features have been 
detected or strongly constrained [2-5]. TRAPPIST-1 b is the 
closest planet to the system’s M dwarf star, and it receives 4 
times as much irradiation as Earth receives from the Sun. This 
relatively large amount of stellar heating suggests that its ther- 
mal emission may be measurable. Here we present photometric 
secondary eclipse observations of the Earth-sized TRAPPIST- 
1 b exoplanet using the F1500W filter of the MIRI instrument 


Springer Nature 2021 TeX template 


2 TRAPPIST-1 b thermal emission 


on JWST. We detect the secondary eclipse in each of five sep- 
arate observations with 8.7-sigma confidence when all data are 
combined. These measurements are most consistent with re- 
radiation of the TRAPPIST-1 star’s incident flux from only 
the dayside hemisphere of the planet. The most straightforward 
interpretation is that there is little or no planetary atmosphere 
redistributing radiation from the host star and also no detectable 
atmospheric absorption from carbon dioxide (CO2) or other 
species. 


The TRAPPIST-1 system has an age of 7.6 + 2.2 Gyr [6] and consists of a 
very cool (Teg = 2566 K), low-mass star (0.09 solar masses) and seven transit- 
ing planets that are 0.75 — 1.1 Earth radii in size with masses 0.3 — 1.4 times 
that of Earth [1, 7]. The small size of its star and its nearby location (12 par- 
secs from Earth) makes its planets amenable to observational characterization 
studies. 

Late M-dwarf stars like TRAPPIST-1 [8] are observed to have frequent 
X-ray and ultraviolet flares [9, 10]. Furthermore, these low-mass stars evolve 
through an extended high-luminosity pre-main-sequence phase lasting 1 Gyr or 
more [11], causing extreme water loss on planets that would later become tem- 
perate [12, 13]. The TRAPPIST-1 planets are also expected to be tidally locked 
so that their rotations are synchronized to their orbital periods [1]. These phe- 
nomena will impact the atmospheres and climates of the TRAPPIST-1 planets, 
causing them to be significantly different from the terrestrial planets in our 
Solar System [14, 15]. High-energy stellar flares and coronal mass ejections 
from M dwarfs can also dissociate molecules and drive atmospheric escape 
[16, 17]. This can be mitigated by planetary magnetic fields if the ejected stel- 
lar plasma heats planetary interiors and leads to volcanism and outgassing 
[18]. Tidal synchronization of planetary rotations to orbits will also cause some 
internal heating that could lead to significant atmospheric outgassing [19]. Fur- 
thermore, tidal synchronization will cause the TRAPPIST-1 planets to have 
warm permanent day sides that always face the star and colder night sides 
that always radiate to free space. 

A number of studies have modeled possible TRAPPIST-1 planetary 
atmospheres, usually with outgassed secondary compositions that may have 
undergone further processing via desiccation or reduction. These atmospheres 
are often dominated by H20, O2, N2, or CO% and have higher mean molecular 
weights than primordial hydrogen-dominated ones [15, 19-21]. TRAPPIST- 
1 planetary transmission spectroscopy observations to date have been precise 
enough to detect only moderately high concentrations of atomic or molecular 
constituents in clear hydrogen-dominated atmospheres [2-5]. These existing 
observations are insensitive to detecting spectroscopic absorptions in cloudy 
or high mean-molecular weight atmospheres of planets as small and cool as 
ones in the TRAPPIST-1 system. 
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The recently-commissioned JWST should be more sensitive to spectral fea- 
tures and thermal signatures of TRAPPIST-1 planetary atmospheres due to 
its improved aperture, wavelength range, spectral resolving power, and stabil- 
ity over that of previous observatories [19, 22, 23]. The inner TRAPPIST-1 
planets are warm enough to potentially detect their thermal emission with mid- 
infrared secondary eclipse observations [24]. TRAPPIST-1 b, the innermost 
planet, has the highest expected equilibrium temperature Teq = 391-400K 
[1, 25]. Teq is computed for absorbing all incident stellar radiation on its day- 
side and re-radiating it isotropically over 47 steradians. Measuring the actual 
dayside temperature would constrain the existence, composition, and circula- 
tion of its atmosphere [26, 27]. Spitzer Space Telescope observations did not 
detect the secondary eclipse of TRAPPIST-1 b but did constrain its dayside 
temperature to below 611 K (3 o upper limit [28]). 

We now report new JWST observations that detect the TRAPPIST-1 b 
secondary eclipse, measure the planet’s dayside temperature, and constrain 
the properties of any atmosphere it may have. Our program consisted of pho- 
tometric observations of five secondary eclipses of TRAPPIST-1 b using the 
F1500W filter of the JWST Mid-Infrared Instrument (MIRI) [29]. This filter 
transmits wavelengths over a half-maximum bandpass of 13.5 — 16.6 microm- 
eters, nearly ideal for maximizing the signal-to-noise ratio of the secondary 
eclipse for the planet’s expected T = 400 K dayside temperature. Each obser- 
vation consisted of 377 integrations of 14 FASTIR groups in full-frame mode, 
with a duration of 4.36 h. Secondary eclipses were observed on 2022 November 
8, 12, 20, 24, and December 3 UTC. The dates and times of the observations 
were computed from a dynamical model of the TRAPPIST-1 planetary system 
[7] and chosen so that there was no contamination by transits or secondary 
eclipses of other planets in the system. 

We reduced the data using the Eureka! pipeline which is optimized to min- 
imize noise in JWST time-series data [30]. Secondary eclipse light curves were 
generated and fit from each observation, and the light curves from all obser- 
vations were also jointly fit (Methods). The secondary eclipse was detected in 
each observation, and Figure 1 shows the combined light curve and fit for all 
observations. 

Next, we convert the measured eclipse depth, F,/F, = 861 + 99 ppm, to 
the measured dayside planetary flux (Fp) and to its brightness temperature 
within the JWST MIRI F1500W photometric bandpass. We used calibrated 
JWST data products to compute a mean stellar flux of 2.589 + 0.078 mJy 
(Methods). Therefore, TRAPPIST-1 b has a dayside flux Fp = 2.229 + 0.263 
uJy. This corresponds to a blackbody brightness temperature Tg = 50313° 
K, computed from Planck’s law. We calculated this Tg value for the F1500W 
filter’s effective wavelength of A = 14.79 um. This wavelength was computed by 
weighting JWST’s wavelength-dependent throughput in this filter by the flux 
of a PHOENIX stellar model [31] interpolated to Tag = 2566 K to represent 
the stellar spectrum. 
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Fig. 1: Combined TRAPPIST-1 b MIRI F1500W secondary eclipse light curve. a. 
The phase-folded light curves from all of our five observations after the removal of 
systematic noise based on Joint Fit #1 (see Methods). Overplotted are data binned 
to a cadence of 9.7 minutes to more clearly visualize the detection of the eclipse, and 
the red line shows the fitted astrophysical model. The model fit has eclipse depth 
Fp/F, = 861 + 99 ppm. b. The residuals after fitting the observations. Error bars 
show lo uncertainties in both panels. 


The apparent observed dayside temperature of a planet in a non-eccentric 
orbit with zero heat redistribution to its night side is Ty = [2/3(1 — A)]'/4T 
where Ty = Tog (R,/a)'/? for semi-major axis a, bond albedo A, stellar effective 
temperature Tig, and stellar radius R, [32]. The planet’s substellar point has 
temperature To, and all other locations have lower temperatures and contribute 
less relative area to the observed planet flux in a secondary eclipse observation. 
For the adopted parameters of TRAPPIST-1 b [7] and A = 0, we compute 
Ta = 508K, very close to our measured Tg = 503738 K. 

Any atmospheric circulation that redistributes stellar heating or absorp- 
tion by atmospheric constituents within the F1500W filter’s spectral bandpass 
will impact the planet’s observed secondary eclipse depth and brightness tem- 
perature. Model TRAPPIST-1 b planetary surfaces composed of anorhite, 
basalt, enstatie, feldspar olivine, pyroxene, quartz, or saponite minerals with 
no gaseous atmospheres are all expected to have secondary eclipse depths 
and brightness temperatures similar to blackbodies in the MIRI F1500W 
bandpass [23]. These and other materials could potentially be identified with 
low-resolution spectroscopic observations at other wavelengths. Furthermore, 
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atmospheres dominated by CO or O> with some CO, (outgassed or des- 
iccated) with surface pressures as low as 10 bar should have significant 
absorptions in the F1500W filter bandpass, reducing the observed brightness 
temperature to 300 K or less when accompanied by efficient heat redistribution 
over the planet [20, 23]. Table 1 shows that such low brightness temperatures 
are inconsistent with our observations. Furthermore, our measured secondary 
eclipse depth is also inconsistent with the Teq = 400 K expected for uniform 
stellar heat redistribution and isotropic re-radiation even with no absorbing 
atmospheric species within the spectral bandpass of our observations. 

The heat redistribution of TRAPPIST-1 b atmospheres have been modeled 
in terms of their surface pressures P and effective gray optical depths 7 [25]. 
Table 1 shows that our observations are consistent with a thin atmosphere with 
little heat distribution; the flux for the surface pressure P = 0.1 bar and gray 
optical depth 7 = 0.1 model [25] is well within 1 o of our measurement. How- 
ever, the atmosphere model with 10 times these values recirculates too much 
heat and has a dayside temperature that is too low to be clearly consistent 
with our observations. 

A coupled atmosphere-interior model that includes heating from tidal 
effects, interior processes, and instellation predicts isotropic surface tempera- 
tures for different TRAPPIST-1 b scenarios at its present age [19]. The median 
isothermal surface temperature of the 50 hottest and lowest albedo PACMAN 
models with negligible final atmospheres (P < 0.1 bar) was 446 K [19] (and J. 
Krissansen-Totton private communication). We converted this mean isotropic 
surface temperature to an apparent observed dayside temperature Tgay = 534 
K using the geometric and thermal conventions for zero redistribution of stel- 
lar heating to the night side as was done in defining Ty above [32]. This should 
be regarded as the maximum Tyay possible from this particular internal heat- 
ing model given that it is derived from only the hottest 10% of model runs 
that resulted in negligible atmospheres. Our observations are consistent with 
this prediction (see Table 1). 


Table 1: TRAPPIST-1 b Observations and Model Comparison 


Atmospheric Model Ref. Psurg(bar) Predicted Tg! (K) Difference (a)? 
96.5% CO2 [23] 10, 92 290 -6.7 
95% O2 + 0.5 bar CO2 [23] 10, 100 303 -6.4 
Isotropic Teq [1] N/A 400 -3.6 
Gray 7 = 1.0 [25] 1 445 -2.1 
Gray r = 0.1 [25] 0.1 495 -0.3 
0 redistribution [32] 0 508 0.2 
0 redistrib. + internal [19] <0.1 534 1.2 


1. Tg in MIRI F1500W filter, when applicable. 
2. Difference between predicted and observed A = 15 um flux in units of measured ø. 


Our observations are most consistent with near-zero bond albedo and little- 
to-no heat redistribution from the dayside to the nightside of TRAPPIST- 
1 b. The most likely explanation for this finding is that the planet absorbs 
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Fig. 2: TRAPPIST-1 b F1500W measured flux and spectral models. The blackbody 
curves represent the measured Tp = 503 K dayside temperature, the 508 K apparent 
dayside temperature predicted for zero heat redistribution and no internal heating, 
and Teq = 400 K temperature for isotropic redistribution of stellar heating. The 
flux expected in the upcoming observations in the MIRI F1280W filter is also shown 
assuming the planet emits like a Tg = 503 K blackbody. The widths of the F1280W 
and F1500W markers represent their transmission half-amplitude bandpasses, and 
the vertical error bar of the F1500W point represents its lo uncertainty. Emergent 
spectra from 93 bar CO% and 10 bar outgassed O with 0.5 bar CO2 atmospheres 
are also plotted [20, 33]. 


nearly all of the incident stellar flux and does not have a high-pressure or 
optically-thick atmosphere. Figure 2 shows the measured TRAPPIST-1 b flux 
in the F1500W filter along with model spectra for blackbodies, CO2, and O2 
+ CO% atmospheres [20, 33]. Our data are clearly incompatible with the Teq 
= 400K, CO2, or O2 + CO2 models. Disfavoring a substantial, high mean 
molecular weight secondary atmosphere suggests that the null results from the 
previous transmission spectroscopy observations [2—4] are likely caused by little 
or no atmosphere of any kind. The absence of a substantial atmosphere is also 
consistent with a coupled atmosphere-interior model that predicts complete 
atmospheric erosion for TRAPPIST-1 b over its lifetime in approximately half 
of the model cases run [19]. 

Future observations have the potential to further constrain the heat redis- 
tribution and possible atmospheres of TRAPPIST-1 b. Spitzer 3.6 um and 4.5 
um band fluxes were once interpreted as possibly due to CO%ə in emission at 
the 2c level [34], but this was not seen in additional observations [28]. JWST 
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secondary eclipse measurements in the MIRI F1280W filter (JWST program 
1279) may better constrain the dayside brightness temperature and determine 
how well the planet’s emission does or does not follow a blackbody emission 
or atmospheric absorption spectrum (see Figure 2). Thermal emission phase 
curve observations that capture most of the planet’s orbit (P = 36 hours) could 
also measure the contrast between its day and night sides. This would allow 
measuring heat transfer and could further determine whether there is a tenu- 
ous atmosphere, as was done with the Spitzer Space Telescope for the much 
hotter rocky planet LHS 3844 b [35]. Secondary eclipse or phase curve ther- 
mal emission observations of additional similar planets including others in the 
TRAPPIST-1 system (e.g., JWST program 2304) could also reveal system- 
atic properties of M dwarf terrestrial planets and how they differ from Solar 
System ones. 
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Data and analysis overview 


These JWST MIRI F1500W time-series observations feature significant back- 
ground signals, dominated by thermal emission from the JWST primary mirror 
and sunshade within the F1500W filter bandpass [36, 37]. At the start of each 
of the five exposures, we measured a mean background of 19,837—22,779 elec- 
trons per pixel in an annulus of radius 12-30 pixels per 38.9 s integration. These 
values assume a gain of 3.1 electrons per digital number (hereafter DN). This 
background level corresponds to approximately 45% of the total flux within 
a five-pixel radius photometric aperture centered on the TRAPPIST-1 star. 
Therefore, we note that it is important to subtract the background flux accu- 
rately to maximize the quality of these and perhaps other long-wavelength 
JWST data. 

We now describe how we analyzed the raw image data to measure the 
secondary eclipse parameters of TRAPPIST-1 b in two different ways. The 
descriptions below refer to calibration pipelines and other software whose code 
and citations appear in the ‘Code availability’ section. 


Eureka! Analyses: Images to light curves 


We reduced our observations using version 0.9 of the Eureka! data analysis 
pipeline [30], with jwst version 1.8.3 [38], CRDS version 11.16.16 [39], and 
using the “jwst_1019.pmap” CRDS context. We started from the _uncal.fits 
files and ran Stages 1-4 of Eureka!. Stage 1 was first run using the default jwst 
settings with the exception of increasing the jump step’s cosmic ray rejection 
threshold to 10 sigma to avoid an excessively large number of false positive 
cosmic ray detections. The only change to the Stage 2 defaults were in turning 
off the photom step except when computing the calibrated stellar flux; turning 
off the photom step makes it easier to understand how close our reduction is 
to the photon noise limit. 

In Eureka!’s Stage 3, we performed aperture photometry using source radii 
ranging from 4-12 pixels. In some of our observations, we found that using a 
large aperture (e.g. a 9 pixel radius which encircles the first Airy ring) resulted 
in a smaller eclipse depth as well as noisier residuals; we attributed this to 
imperfect background subtraction across the large source aperture resulting in 
a dilution of our eclipse depth. Ultimately, we found that an aperture radius 
of 5 pixels (which encircles the Airy disk) offered the lowest scatter in our 
residuals by maximizing the encircled starlight and minimizing the encircled 
background flux. We also considered a range of annuli to determine the back- 
ground in each frame which was important as the background flux varies 
significantly over the duration of the observations. We found that the choice 
of background annulus had little impact on the resulting eclipse depth, but 
an annulus that started near the end of the first Airy ring and included many 
background pixels helped to reduce the scatter in our residuals; ultimately we 
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chose an annulus spanning 12-30 pixels from the centroid. We also masked 
the values marked as DO NOT USE in the DQ array, 50 clipped background 
flux outliers and linearly interpolated bad pixels. We computed the centroid 
for each integration using a 2D Gaussian fit with an initial value for the cen- 
troid of x=698, y=516 and using only the 20x20 grid of pixels surrounding 
this first value. In Eureka!’s Stage 4, we sigma clipped 100 outliers compared 
to the local median flux computed using a boxcar filter with a width of 10 
integrations. 

We also performed an alternative data reduction that was identical to the 
one described thus far except we set the ramp_fitting_weighting to uniform 
in Stage 1. All groups in each integration of the data are limited by photon 
noise (background and star signal) and not read noise, so it is logical to weight 
them uniformly instead of using the default optimal value which is optimized 
for faint sources and low backgrounds. 


Calibrated Stellar Flux 


We measured the stellar flux of TRAPPIST-1 from _calints.fits data prod- 
ucts that we downloaded from the JWST archive. These had undergone 
standard processing with calibration version 1.8.2 and data processing version 
2022_4a. After discarding the first 40 integrations of each observation, we com- 
puted the flux in a 15 pixel radius aperture and subtracted the background in 
a 16-35 pixel annulus centered on the star. We applied corrections for starlight 
outside of the measurement aperture and inside the background annulus (~4% 
each; [40, 41]) and applied the pixel scale of 2.844 x 10713 sr/pixel. The cali- 
brated flux varied by < 0.4% (o = 0.004 mJy) between the five observations, 
and there is an additional absolute calibration uncertainty ~3% (0.078 mJy; 
K. Gordon / Space Telescope Science Institute private communication). We 
add the systematic and observational uncertainties in quadrature to get an 
observed stellar flux of 2.589+0.078 mJy. This is 13% lower than than the flux 
expected from a Phoenix / BT-Settl model produced by the pysynphot pack- 
age for the stellar parameters and normalized to the star’s observed 2MASS J 
= 11.35 mag flux. 


Secondary eclipse depths 


To begin, we first individually fit each of our visits to assess the data quality of 
each visit and assess the evidence for astrophysical variability. We performed 
two independent fits to the observations to ensure the reproducibility of our 
conclusions. Individual and joint fit #1 used the standard-processed Stage 
1-4 data, while individual and joint fits #2 used the alternatively processed 
data with ramp_fitting_weighting=uniform. We describe the method used 
in each fit below. 
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Individual Fits #1 


For the first fit to the individual observations, we used Eureka! to fit a uniform- 
surface eclipse depth model using starry [42], a linear trend in time, and we 
linearly decorrelated against the star’s x and y position as well as the star’s 
x and y PSF width as determined using the Gaussian centroiding method in 
Eureka!’s Stage 3. We also inflated the estimated white noise level during the 
fit to account for noise above the photon limit. The first ~100 integrations 
showed a moderate ramp shape that was likely caused by detector settling, and 
we decided to remove these data rather than try to fit them with an exponential 
decay model. To fit the observations, we used PyMC3’s No U-Turns Sampler 
[43] and ran two chains with a target acceptance rate of 0.95 with each chain 
taking 5000 tuning steps and then 2000 posterior samples. The Gelman-Rubin 
statistic [44] was used to ensure the chain had converged. For astrophysical 
priors, we took R,/R, = 0.08590, R, = 0.1192 Re (to account for light travel 
time), F,/F, = 0+ 0.01, P = 1.5108794 days, to = 59890.5150313 + 0.005 
(BMJD_TDB), i = 89.728°, a/R, = 20.843, and e = 0 [7]. From these inde- 
pendent fits, only the eclipse from November 20th showed evidence for a small 
amount of unmodelled red noise in its Allan variance plot [45]; this unmod- 
elled red noise is likely a slight downward ramp at the start of this observation 
(see Extended Data Figure 1). 

Note that our to prior comes from numerical predictions which account 
for planet-planet interactions which significantly perturb the planets’ orbits 
from static, circular orbits [7]. To ensure that any inaccuracies in those pre- 
dictive models do not bias our eclipse depth, we inflated the uncertainty on to 
compared to the published uncertainty [7] to search through a larger region 
of parameter space. Reassuringly, all of our individual eclipse timings were 
within error of each other and within 2 ø of the predicted timings [7] assum- 
ing zero eccentricity. Each observation constrained the eclipse timing to within 
~3 minutes, or ~8% of the eclipse duration. The mid-eclipse times for the five 
observations are given in Extended Data Table 1. 

Our five independently fit observations gave eclipse depths of 790 + 210 
ppm, 510+210 ppm, 950+170 ppm, 820+ 220 ppm, and 829+ 200 ppm. Each 
of these observations is consistent with an average eclipse depth of 795 ppm, 
and we find no evidence for variability in the eclipse depth of TRAPPIST-1b. 


Individual Fits #2 


For the second fit to the individual observations, we used the Fortran code 
trafit which is an updated version of the adaptive Markov-Chain Monte 
Carlo (MCMC) code described in [46-48]. It uses the eclipse model of [49] 
as a photometric time-series, multiplied by a baseline model to represent the 
other astrophysical and instrumental systematics that could produce photo- 
metric variations. First, we fit all visits individually. We tested a large range 
of baseline models to account for different types of external sources of flux 
variations/modulations (instrumental and stellar effects). This includes poly- 
nomials of variable orders in: time, background, PSF position on the detector 
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(x,y) and PSF width (in x and y). Once the baseline was chosen, we ran a pre- 
liminary analysis with one Markov chain of 50,000 steps to evaluate the need 
for re-scaling the photometric errors through the consideration of a poten- 
tial under- or over-estimation of the white noise of each measurement and 
the presence of time-correlated (red) noise in the light curve. After re-scaling 
the photometric errors we ran two Markov chains of 100,000 steps each to 
sample the probability density functions of the parameters of the model and 
the system’s physical parameters, and assessed the convergence of the MCMC 
analysis with the Gelman & Rubin statistical test [44]. Our jump parameters 
were the following : the mass of the star M, with a normal prior distribu- 
tion N(0.0980, 0.00237) Mo; the radius R, with a normal prior distribution 
N (0.1192, 0.00337) Ro; the effective temperature Tepp with a normal prior dis- 
tribution M (2520, 392) K; and the metallicity [Fe/H] of the star with a normal 
prior distribution \(0.0535, 0.08807); the impact parameter of the planet b 
with a normal prior distribution \(0.25,0.117); and Tọ with normal priors 
and inflated errors similarly to Fit #1; and F,,/F, with a broad uniform prior. 
All normal priors were taken from [28] except for the transit timing that were 
derived from the dynamical model predictions by [7]. All of our individual 
eclipse timings were also within 2 ø of this dynamical model’s predicted timings 
assuming zero eccentricity. 

For the five individual fits the Gelman-Rubin statistic was lower than 1.003 
for every jump parameter measured across the two chains, indicating good con- 
vergence. The eclipse depths that we measured for visit 1 to 5 are: 9947183 ppm, 


691; 168ppm, 7987382ppm, 8217225 ppm, 7367 2gppm. 


Combining the multiple eclipse measurements 


Finally, we simultaneously fit all five of the eclipse observations to allow for 
better constraints on the timing of eclipse; this is important as any error in 
the mid-eclipse time will result in a bias toward smaller eclipse depths as some 
of the eclipse signal will included in the baseline and vice versa. 


Joint Fit #1 

For the first joint fit to all observations, we followed the same method as 
the Individual Fits #1 with the exception of fitting the period using a Gaus- 
sian prior of P = 1.5108794 + 0.000006 [7] and assuming a shared F,/F,, P, 
and to between all visits; all other fitted variables were allowed to vary freely 
between the different visits. We again fit the observations using PyMC3’s No 
U-Turns Sampler [43] and ran two chains with a target acceptance rate of 0.90 
with each chain taking 6000 tuning steps and then 3000 posterior samples. 
Again, the Gelman-Rubin statistic [44] was used to ensure the chain had con- 
verged. From this joint fit, again only the eclipse from November 20th showed 
evidence for a small amount of unmodelled red noise in its Allan variance 
plot [45]. This shared fit constrained the eclipse time to within 36 seconds 
(to = 2459891.01487704778 + 0.00042; or ~1% of the eclipse duration) and 
is in close agreement with the numerical predictions [7]. Our fitted orbital 
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period (P = 1.5108793 + 0.0000062 days) also remained very consistent with 
the prior value [7]. Overall, we find that the numerical, planet-planet inter- 
action model [7] based on observations taken three or more years before our 
JWST observations was able to predict the timing of TRAPPIST-1b’s eclipses 
with excellent accuracy and precision. We also find that the eclipse duration 
is consistent with the planet’s published transit duration, but our limited data 
allow constraining this time to only ~10%. 

The individual raw and de-trended light curves as well as the joint model 
fits to our phase-folded observations are shown in Extended Data Figure 1. 
With this joint fit, we measure an average eclipse depth of 899+91 ppm. While 
this average eclipse depth is significantly higher than that from the individual 
fits, this seems to be the result of a bias towards smaller eclipse depths in the 
individual fits as a result of higher uncertainty on the eclipse timing. Indeed, 
allowing the eclipse depths to vary between visits while only requiring the time 
of eclipse to remain consistent between visits gives an average eclipse depth 
entirely consistent with our joint fit. The phase-folded light curves from all 
observations after systematic noise removal are shown with the joint eclipse 
model fit in Fig. 1 (Main). 


Joint Fit #2 

For the second joint fit to all observations, we used trafit with the baseline 
models derived from our individual fit for each light curve (see Individual Fits 
#2). Again we did a preliminary run of one chain of 50,000 steps to estimate 
the correction factors that we then apply to the photometric error bars and 
then a second run with two chains of 100,000 steps. The jump parameters were 
the same as for the individual fit #2 expect from the fact that we fixed To and 
allowed for transit timing variations to happen for each visit (each transit TTV 
has an unconstrained uniform prior). We used the Gelman-Rubin statistic to 
asses the convergence of the fit. We measure an eclipse depth of 82373 ppm 
from this joint fit. 


Adopted Eclipse Depth and Noise 


We compute a single mean secondary eclipse depth F/F, = 861 ppm from the 
two separate jointly-fit values of 899 and 823 ppm. We adopt the uncertainty 
99 ppm = (ø? + 03)!/? where cı = 91 ppm is the larger of the two jointly-fit 
random uncertainties, and 02 = 38 ppm is an estimate of the systematic error 
equal to the difference between either jointly-fit depth and their mean. The 
random noise value of 91 ppm is ~1.3 times the photon noise we measure from 
the data and also estimate from an in-flight JWST MIRI performance model. 
This excess noise factor is similar to one measured for the MIRI low resolution 
spectroscopy (LRS) time-series observing mode which uses the same detector 
and signal-chain electronics [50]. 
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Data Availability 


The data used in this paper are associated with JWST GTO program 1177 
(observations 7 — 11) and will be publicly available from the Mikulski Archive 
for Space Telescopes (https://mast.stsci.edu) at the end of their one-year 
exclusive access periods. 


Code Availability 


We used the following codes to process, extract, reduce and analyse the data: 
STScI JWST Calibration pipeline [38], Eureka! [30], Emcee [51], starry [42], 
PyMC3 [43], PySynphot [52], and the standard Python libraries numpy [53], 
astropy [54], and matplotlib [55]. These were incorporated into custom python 
notebooks for data analysis. These notebooks are available on request. The 
notebooks were developed by a NASA employee and are not available publicly 
until approved by NASA. 
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Extended Data 
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Extended Data Figure 1: Light curves of individual observations and the Joint 
Fit #1 models. a. The raw light curves from each of our five visits normalized by 
their median value are shown in color, while the fitted model from Joint Fit #1 
including systematic noise is shown using a black line. The date in UT of each visit is 
indicated by the y-axis labels. Each visit did not start at the same orbital phase, so 
the apparent movement in the eclipse time is simply caused by when the observations 
began. b. The same data and model from each visit after the removal of systematic 
noise. Overplotted are data binned at a cadence of 9.7 minutes (14 integrations) to 
more clearly visualize the detection of the eclipse in each visit. All error bars show 
lo uncertainties in both panels. 
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Extended Data Table 1: Mid-eclipse Times 


Observation UT date (2022) mid-eclipse time! (BJD TDB) 


1 8 November 2459891.771179-bo16 
2 12 November  2459896.303848-0819 
3 20 November 2459903.855919:9011 
4 24 November 2459908.3893%0-0012 
5 3 December —2459917.4563 40-0090 


1. Stated uncertainties are 1 ø. 


